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Abstract 

We construct a one-particle TOA operator T canonically conjugate 
with the Hamiltonian describing a free, charged, spin-0, relativistic par¬ 
ticle in one spatial dimension and show that it is maximally symmetric. 
We solve for its eigenfunctions and show that they form a complete and 
non-orthogonal set. Plotting the time evolution of their corresponding 
probability densities, it implies that the eigenfunctions become more lo¬ 
calized at the origin at the time equal to their eigenvalues. That is, a 
particle being described by an eigenfunction of T is in a state of definite 
arrival time at the origin and at the corresponding eigenvalue. We also 
calculate the TOA probability distribution of a particle in some initial 
state. 


1 Introduction 

The notion of time is one of the problems obstructing the marriage of Einstein’s 
General Theory of Relativity and Standard Quantum Mechanics into one frame¬ 
work of Quantum Gravity [I]. The two theories have a mutually incompatible 
treatment of time. For general relativity, time has a dynamic and intrinsic role 
in the evolution of the system being studied. For quantum mechanics, however, 
time is merely an extrinsic parameter marking the evolution of the system. The 
system does not affect it nor does it affect the system. This pessimistic view on 
time was prominent in the earlier days of quantum mechanics. The prevalent 
formulation was the von Neumann (standard) formulation of quantum mechan¬ 
ics which implied that time is not a dynamic physical entity to be measured 
- time, according to that view, is not an observable. Observables in quantum 
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mechanics are represented by mathematical objects called operators and in the 
standard formulation, these operators should be self-adjoint. The demotion of 
time into just being a parameter stems from Pauli’s theorem which says that 
there can be no self-adjoint time operator canonically conjugate with its corre¬ 
sponding semi-bounded system Hamiltonian [2]. This pessimism towards time 
was apparent through the trilogy of papers by Allcock ElllIS] which suggests 
that one cannot find a distribution for the time it takes for a free particle to 
arrive at a certain point - the particle’s Time of Arrival (TOA). 

This pessimism seems misguided however, since we can, for example, mea¬ 
sure the TOA of a particle. As an answer to this, studies to promote the role of 
time has been made in recent years. In the framework of Positive Operator Val¬ 
ued Measures (POVMs), the standard formulation of quantum mechanics can be 
extended by saying that quantum observables are not necessarily self-adjoint op¬ 
erators but may be non-self-adjoint, maximally-symmetric operators which are 
the first operator moments of POVMs. This framework lets one bypass Pauli’s 
theorem by allowing one to consider non-self-adjoint, maximally-symmetric time 
operators as quantum observables. Moreover in [^, it was rigorously shown that 
the assumption of the existence of a bounded, self-adjoint time operator canon¬ 
ically conjugate to a semi-bounded Hamiltonian is consistent. It was seen that, 
in the steps followed by Pauli, there were some implicit assumptions made and 
that these assumptions were inconsistent. It is then the case that Pauli’s the¬ 
orem simply does not hold in the standard formulation of quantum mechanics. 
This then opens up an avenue to still consider time as an observable in standard 
quantum mechanics. In line with this optimistic view, several studies on TOA 
operators [H |8l El HOl HH 1121 UHl IHl IH IHl Ull IHl US] were made. Specifi¬ 
cally in [12], it was seen that a time operator, conjugate with its corresponding 
Hamiltonian, can be written as a sum of the Bender-Dunne basis operators Tm,n 
jiniEl]. As an example, the time of arrival for a free non-relativistic particle in 
one spatial dimension would then be given by — [Z11I2], where m is the 
mass of the particle. Some applications of these TOA operators would be a test 
to the possibly inconsistent assumptions of neutron time-of-flight spectroscopy 
experiments [23] (that is, simultaneously assuming the validity of the classical 
time of arrival and the quantum mechanical broad de Broglie wave packets), 
the delay in photoemission from atoms which is in conflict with the assumption 
of spontaneous photoemission [24] , and the upper limit to a quantum tunneling 
delay time being much shorter than previously predicted [25] . 

In our quest to promote time as an observable, that is, to make time be on 
equal footing as the other physical quantities, the next system to study time in 
would be the quantum mechanics of a relativistic particle [n EH na im |30]. 
In this regime, special relativity becomes relevant and it tells us that space and 
time should be treated on the same ground. If position is to be taken as an 
observable of the particle, it is then natural to promote time (such as the TOA 
at a certain location) also as an observable with a corresponding operator. This 
TOA operator should then be conjugate with the system Hamiltonian which 
may be derived from the equations of motion describing the system. Another 
property of the relativistic quantum particle is its spin. Particles with different 
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spins are described by different equations of motion. Some of the well known 
equations of motion are the Klein-Gordon equation, describing spin-0 particles, 
and the Dirac equation, describing spin-1/2 particles. However, there are cer¬ 
tain problems arising in these equations, problems such as ’negative energies’ 
and ’negative probabilities’. These can be attributed to the existence of anti¬ 
particles. They are just particles with the same mass but opposite charge of their 
corresponding particle pair. The particle anti-particle pairs are spontaneously 
created and annihilated in a relativistic quantum system so that the system 
does not have a fixed number of particles. This should not be too surprising 
since quantum mechanics tells us that energy can have quantum fluctuations 
and special relativity tells us that mass is just another form of energy, then 
it is natural that relativistic quantum systems may have a fluctuating number 
of massive particles. Standard quantum mechanics, however, was not formu¬ 
lated with a changing number of particles in a system. Quantum Field Theory 
(QFT), the theory of quantum fields where space and time are just labels of a 
point on the field, remedies this by saying that particles and anti-particles are 
just excitations of the quantum field which may undergo random quantum fluc¬ 
tuations, thus the creation and annihilation of particles and anti-particles. This 
is one of the strengths of QFT. We can, however, still ask the question: Given 
a relativistic particle with a certain initial state, what is the TOA probability 
distribution of that same particle? If we are to perform this TOA experiment, 
we would not be sure if the particle that arrived at our detector is the same one 
we started with. But surely, we need not the full machinery of QFT to answer 
the quantum mechanics of this relativistic particle. We just need to restrict our 
analysis to one particle [33]. That is, we would only consider the cases where it 
is the same particle that has arrived at the detector of our TOA experiment. 

The purpose of this paper then is to construct a TOA operator for a free, 
relativistic, spin-0, charged particle in one spatial dimension by solving for an 
operator T conjugate with the system Hamiltonian derived from the Klein- 
Gordon equation. We are also explicit in making the constructed TOA operator 
a true one-particle operator. Other papers have also studied time operators in 
a relativistic context. These studies include |5Bj which has also constructed a 
relativistic TOA operator for spin-0 particles which closely mirrors [5] . That is, 
the time evolution of the position operator was inverted to obtain an expression 
for a TOA operator defined with a suitable operator ordering. In [27], however, 
a kinematic time operator was defined for relativistic quantum mechanics but 
the dynamics suggest that the eigenstates of this time operator is unphysical. 
To remedy this, the Bohmian interpretation was invoked. While in [55], a self- 
adjoint time operator was constructed for a free massless relativistic particle in 
terms of the Poincare group generators. Under Lorentz boosts, the transforma¬ 
tion of the operator differs from its expected transformation law from special 
relativity. This suggests that the concept of time associated with this operator 
differs from the Minkowski time coordinate. Another is in [55], where an alter¬ 
native set of coordinates are considered to avoid the false result of having the 
eigenvalues of the velocity operator of a massive relativistic particle be ±c which 
suggest that particles with mass travel at the speed of light. The corresponding 


3 


operator of the ’zeroth’ component of these new coordinates does not commute 
with the Hamiltonian of the system which gives us a time-energy uncertainty 
relation. In |30] . a self-adjoint dynamic time operator for relativistic spin-1/2 
particles was constructed. The resulting commutation relation with the system 
Hamiltonian is analogous with the position-momentum commutation relations. 
Such relativistic time operators would be useful in atomic fountain clocks where 
they lower the temperature of the particles to hinder relativistic effects m- 
These procedures would no longer be needed if one can study the quantum be¬ 
havior of time in a relativistic regime. One can see then that there are various 
approaches in the study of time in relativistic quantum mechanics. The hope is 
that upon further study, these approaches would give us a more refined under¬ 
standing of time. Moreover, promoting time to be on the same footing as other 
physical observables may give us the insight to better understand the problem 
of time in Quantum Gravity [I]. 

In section [51 we first review the one-particle interpretation of the Klein- 
Gordon equation by closely following the discussions from [33) . Here, we review 
how we can separate wavefunctions to describe either a positively charged or 
negatively charged particle, allowing us to have a non-negative probability den¬ 
sity of the charged particle and a notion of true one-particle operators. We 
construct the TOA operator by solving its canonical commutation relation with 
the Hamiltonian of the Klein-Gordon particle in section [31 make it into a true 
one-particle operator (which we will later call as T) in section 21 and show that 
it is symmetric in section |5l We solve for the eigenfunctions of T in section |6l 
show that they form a complete and non-orthogonal set in sections [7] and [HI re¬ 
spectively, and study their dynamical behaviour by plotting the time evolution 
of their associated probability densities in section |9l In section [101 we calculate 
the probability (density) that the particle in some initial state will arrive at the 
origin at time r and, as an example, considered the initial state to be normally 
distributed about some initial position and momentum. Lastly in section 1111 
we conclude. 


2 One-Particle Interpretation of Klein-Gordon 
Particles 


In one spatial dimension, the wavefunction tjj oia, free, relativistic, spin-0 particle 
with mass mo and magnitude of charge e is described by the Klein-Gordon 
Equation (KGE) 


1 d'^ip 9 ^^ niQC^ 

dt"^ dx^ h? 

Along with Eq. ([T]), a 'probability density' p can be constructed 


( 1 ) 


ih 

2moc^ 




( 2 ) 


However, since the KGE is second order in time, both ip and dip/dt can have any 
arbitrary value at some given time. This means that p is not positive definite 
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and hence, not a probability density. This can be remedied by interpreting it 
as a charge density. That is, p = ep. In this interpretation, p can then be 
allowed to be positive or negative as it just measures the difference between the 
number of positively charged and negatively charged particles. One sees then 
that there is another degree of freedom in ?/; corresponding to states describing 
positively and negatively charged particles. To be more concrete, consider the 
plane wave solutions for a particle with definite momentum p which has the 
form '0 oc exjp(ipx/h— iEt/h). Then substituting i/) back into Eq. ([T]), we get 
E = = ±Ep = XEp where, A = ±1 and Ep = ^/p^c^ + m^c^, 

so that for a particle with a definite momentum p, we get two 'energies' giving 
rise to two different plane wave states. We can then characterize these plane 
wave states by the sign of A. Explicitly, the plane wave states are '0 a=±i c>c 
expippx/h — XiEpt/h). Moreover, the corresponding charge density p = ep 
of '0A=+i ('0A=-i) is positive (negative) definite. We can then interpret that 
■0A=+i {ip\=-i) describes a positively (negatively) charged particle. Note that 
E = ±Ep = XEp does not imply that there is a positive and negative energy. 
The sign of A simply characterizes whether the state describes either a positively 
charged or negatively charged particle. 

However, as presented in Eq. o, the KGE does not make this separation 
of solutions into positive and negative states obvious. It is best to make the 
charge degree of freedom more visible by decomposing the second-order-in-time 
KGE (Eq. (Ill) into two first-order-in-time differential equations. Moreover, 
this allows us to define a Hamiltonian operator for the system. The first natural 
step is to let 0 and dtp/dt to be independent functions. Explicitly from [33], we 
let ijj = (p + X and ih^ = moc^{ip — x) so that Eq. dH becomes separated into 
two equations made compact by a two-element column vector = ((/? y)^. 




dt 


\ 02 


ih— = (ct 3 -b ia2) --— - 5-2 ^ + CT 3 TO 0 C 


2mo J dx"^ 


where, is the Hamiltonian operator of the system, and 


CTO = 


1 0 
0 1 


,CTl = 


0 I 
1 0 


,CT2 = 


0 -i 
i 0 


,CT3 = 


I 0 
0 -1 


(3) 


(4) 


Note that ctq is just the identity matrix. We will follow [33] and call Eq. m 
the Schrodinger representation of Eq. dn, or the ^/-representation. The charge 
density can also be written as 


p = ep = e{ip*ip - X*x) = <^3'^ (5) 

Eq. d3|: however, is still a coupled differential equation and we can decouple 
it by diagonalizing In the case we are investigating, the non-interacting 
case, this is possible. Firstly, it is more convenient to work in the momentum 
representation (p-representation) so that we make the replacement —ihd/dx —> 
p in the Hamiltonian T-L^. The functions from this point will also be in the 
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^-representation unless otherwise stated. Then from [311132] ) we introduce 


jj±i _ (wqc^ + Ep)aQ =F (mpc^ - Ep)ai 


(6) 


where 17 ^ is the inverse of U. The KGE (Eq. ([3|)) then takes the form 

= E;p0-3$ = 77$$ (7) 

where $ = 175' = (</)+ is a two-element column vector and 77$ is the 

corresponding Hamiltonian. Following |33j . we call this representation the 
Feshbach-Villars representation or ^-representation for short. Also, if we want 
to transform an arbitrary operator from the '^-representation to the <f>- 
representation, we simply use the transform = 1771^17“^. We sometimes 
indicate $ — p representation to emphasize that we are also in the momentum 
representation. Note that since 77$ is diagonal, Eq. © does not mix </)+ and 
(j)-. That is, the free time evolution of </>+ (</>_) depends only on </)+ {(j)-) itself. 
To be more concrete, if we have $+ = (</)+ 0)^ and $_ = (0 then Eq. ([7|) 
becomes 

Ep(l>+ ( 8 ) 

-Epct>_ (9) 


ih 


ih 


d(t>+ 

dt 

d4>- 

w 


Moreover, we can calculate the charge density (albeit, density in momentum 
space in our present discussion). Following from Eq. ([5|) 

p = e'I'Vs'I' = e<i>V3$ = e (|</>-i-|^ — (10) 

so that the corresponding charge density for $-|_ ($_) is e|(()+p {—e\(j)-\^). That 
is, $+ ($-) represents a positively (negatively) charged particle. One sees then 
more explicitly that had we started with a free positively (negatively) charged 
particle, then after its evolution given by Eq. ([5|) (Eq. ©), we still get the same 
positively (negatively) charged particle. With this, restricting our analysis to 
one particle is more straightforward. 

Since the corresponding charge density of $+ ($-) is positive (negative) 
definite, we can treat \4>±\'^ as a probability density (again, density in momentum 
space in our present discussion). That is, from the total charge of the positively 
charged or negatively charged particle 


/ P±dp = ±e 

J — OO 
>o 

$)j_cr3$±(7p = ±1 

CX3 

/ OO 

\4'±\^dp = 1 

-OO 
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we see that |</>±P is normalized to unity. Also, the second line can be interpreted 
as the normalization of the state <1>±, or < >$= ±1. This implies that 

we have an inner product defined as < >$= ^\a 3 ^ 2 dp for some 

and $2 in the (implied) system Hilbert space of states H 

< oo (11) 

Restricting our analysis to one particle then is much like splitting the Hilbert 
space into the positive and negative states and working in one of the subspaces. 
That is, we have H = 71+ 0 H- where the elements of 77+ {77-), in the $- 
representation, all have the form $+ = (0+ 0)^ (<&_ = (0 4>-)’^) so that if we 
are studying a positively (negatively) charged particle, then we only work on 
the Hilbert (sub)space 77+ {77-). 

We may also consider operators in this Hilbert space 77. In general, an 
operator A may be decomposed as a sum of a diagonal operator and a non¬ 
diagonal operator A=[A\ + {A} where, [A] is the diagonal operator (so-called 
even operator) and {A} is the non-diagonal operator (so-called odd operator) 
[SH [33] . If we are only considering one particle, then true one-particle operators 
should not map any element of 77+ into 77- and vice-versa so that the charge 
sign of the particle is preserved. In the ^-representation, it is easily seen that 
these operators are the diagonal operators so that the true one-particle operators 
are the even operator parts. That is, given a general operator A in the $- 
representation, its even part [A] is the true one particle operator. An example 
of an even operator is the Hamiltonian 77<^ and the momentum operator. The 
position operator, however, is not and we need to take its even part if we wish 
to use it in one particle analysis. 


77 




^£77 


3 Constructing 7^ 

In this section, we attempt to construct a Time of Arrival operator canonically 
conjugate with the system Hamiltonian with the correct non-relativistic limit. 
For simplicity, we begin with the canonical commutation relation 


77-+,'%! 


dt = 


( 12 ) 


where, 77^ is still the Hamiltonian in Eq. ((Sj) in representation-less form. It is 
explicitly given by 


77 ^ = 


1 1 
-1 -1 


2mo 


1 0 
0 -1 


rufjcri 


= {as + 10 - 2 )-^-- 1 - asuioc'^, 

zmn 


( 13 ) 
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the CTj’s are given by Eq. o, and tp is a two-element column vector. We intend 
to solve for the operator by letting 



(14) 


where the Am,nS are unknown but constant 2x2 matrices, and the Tm,nS are 
explicitly given by 




(15) 


which can be extended for either negative m or n [20l [2T]. The Tm.n’s form a 
complete linearly independent set and are just the Weyl-ordered quantizations 
of the one dimensional monomials of the classical position and momenta 
In a recent paper |22| . it was shown that the T-m,nS for positive m and n are 
densely-defined operators in the Hilbert space £^(R) and are thus, meaningful 
quantum mechanical operators. We use the first equality in Eq dTSl) as our 
expression for the Tm,n’s and to simplify our calculations, we let the upper 
limit of the sum over k be positive infinity since the binomial coefficient 
vanishes for k > n. Since the Uj’s are linearly independent and form a complete 
set, we can write any 2x2 matrix in terms of them. Specifically, we can write 


3 



(16) 


where the aj^’^’s are just unknown scalars. 
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Substituting Eq (HI cinci Eq (|16p into tti6 left Imnci side of Eq (jl2p , we Eclvg 




^ = E 


(^3 “1“ ^^ 2 ) 1 

2uIq 


^AYn,nTm,T 




m,nP ~Tm,n Am,n {^3 H“ ^ 


= E 


+moc^ [a^^A 

m,n ] fm,n) 

( 0-3 + ia 2 )A 

m, 




2mo 

Am,n(^^3 A io' 2 ') 

2mo 


Tm+2,n 1 ^ Ti{tI ^^TjYi yi—2 


;i2 . 

Tm+ 2 ,n + irmTm+l,n-l -~ l)Tm,ra-2 


+moC^ [ 0 - 3 , ^m.™] ^m,n^ 


= E 




<73 + J(T2, 7-- Am- 2 ,n ~ 

ZrriQ 


h^{n + 2)(n + 1) 

8mo 


A. 


m,n+2 


tfi(n + 1) 


2mn 


^m —l,n+l} “t“ UIqC [(T3, ^m,n] ) 


= E 


1 m —2,n ^ “t “b 1) m,n+2 \ r)/ , ■ \ 

- 0t^ -- a^’ 2((T3 + lcr2) 


2?71o 


8too 


b m—2,n ^ “b “b 1) m,n+2 

“Q!o QI 9 


2too 


8mo 


b m—2,n h {jl -\- 2)(u + 1) 7 n n+2 

-a,---ao 


2too 


8mo 


2icri 

2(71 


ih{n + 1) 
2mo 


( m —l,n+lo/ I • \ I m—l,n+ln- i m —l,n+lt-\ \ 

ao 2(f73 + 20 - 2 ) + 0(2 2zf7o + 0^3 2f7o 1 


+moc^ (a™’"2icr2 - Q!™’"'2zcri)) T'm.n'I' 


(17) 


where the fourth equality is obtained by shifting indices and some rearrangement 
of terms, while the last equality is the right hand side of Eq (fT^ . Collecting 
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the coefficients of the cr^ ’s, we arrive at four equations for the a™’"’s 


ih{n + 1) 


2m. 


0 


-2 


( m —l,n+lo- I m —l,n+lo\ -tr r 

0^2 2?. 2 1 — ^^(^?7i,0^n,0 

h^{n + 2)(n + 1) 


m—2,n I • m—2,n 

i;q + lao 


X a 


2mo 8mo 

( m,n+2 I • m,n+2 


— 2imQC^a^’^ = 0 


n - ( ^ m—2,n ^ “1“ “1“ _ m,n+2 

I _ Ci 1 _ Ct 1 


2mo 

-2i 


8mo 

ih(^Tl 1 ) ^_|_x . 2 m,n n 

-ttf, + Zim^c an =0 

2mo 


(-) I 1 m— 2 ,n ^ “f “f 1) m,n +2 

2 T-an---a. 


2mn 


-2 


ih{n + 1) 
2mn 


8too 

m —l,n+l 
^0 


= 0 


(18) 


(19) 


( 20 ) 


( 21 ) 


Combining Eqs (uni and m yields a™’” = 0 for any m, n which in turn 
gives = 0 for any m from Eq m- Simplifying Eq m, we get 

n (a™’" + = —moSm,-iSn,i then plugging it in Eq (IT^l) : 


2mo 


fi^n(n + l)mo6m-iSn,-i 
8too 


2 m.n n 

imoc na2 = U 


which gives = (2jmoc^n)“^(5m,i(5n,i and then using Eq (fT51) a™’"^° = 

— {2moc^n)~^Sm,iSn,i — 'moSm,-i^n,i/n for any m. Note that we have some 
arbitrary constants which may not necessarily vanish i.e. a^Q 2 3 0 for any 

m. For simplicity however, we wish to take 7$ as the minimal solution of 
Eq (112[1 . That is, we set as many a^- ’ 's as possible to vanish m so that 
the only non-vanishing constants are explicitly given by a^’^ = {2imoc^)~^, 
ag’^ = —(2 tooc^)~^: and a^^’^ = —mo. 

The minimal solution of Eq (1121) is then 


-1.1 


%, = (^al’^a2 + ag’V3 j Ti^ -(- ag V3r_i4 
(0-3 -I- i(T2)fi^i - moaoT-1,1 


2moc^ 

1 

2moc^ 


pq + qp 


mo 


0 ^ ^ P ^q + qp 


( 22 ) 


so that in the Schrodinger-momentum representation (tk—p representation), we 
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have the action 


(-1 -i) I + Vp 

-“"(J -i)?(;!*<')+ !(;*<>’> 


(23) 


4 Calculation of the even part of 7$ 


We can calculate the $ — p representation of Eq (1^51) by using the transform 
7i$ = for some trial row vector $ where, U and U~^ are given by 

Eq. ([5]). 

Consider first Ti iU~^^ in the momentum representation 


=I (4 c"*)+1 

=1 Kl"'") 




= + ihp ( -u 


p-i PC^cri 

2El 


$ 


= [/- 


2„2 




ihp‘‘c 


■CTi^h 


and similarly, T_i iC7 


(24) 




C/“4 


2 \p \dp J p 


op 


+ ("1$ 

op J p op \p 


= + (-u 


p 


r-i PC 0-1 

2E2 


$ 


= U- 




(25) 
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So that we can calculate: 


1 


2moc^ 


U+ ia2)Ti^iU — m^lJazT-i^iU 


^ 17(CT3 + ifT2)C/-i 


2moc^ 

— moUa3U~^ 
1 


2EI 


T_i,i$- 


2EI 


CTi$ 


2moc^ 




f ' 

— TOn - 


moc 
mo I —f;—CT s 


= 0-3 


V 




2£'2 


2En 


En 


P 

2Er 


T_ 


1,1 


+ <72 




ip f _^ 

2Ep 2Ep3 


$ 

$ 


(26) 


The even part of 7$, denoted by [7$], is the diagonal part of 7$. [7$] 

represents the operator part that does not mix the positive and negative states 
and is therefore a true one-particle operator. Explicitly in $ — p representation, 
its action is given by: 


([r$] $) 


2Ep 
E-n 




\-l 

Vo 

-ij 

' 2 V 


2^2 


777-0 C' 
2Ej3 


2c2 







) Vo 

-ij 

' 2 V 


19 d 


$(p) 


(27) 


Being a one-particle operator, we can interpret 



= T as two indepen¬ 


dent operators, one for positively charged states and another for negatively 
charged states. Considering the operator for positive states and taking the non- 
relativistic limit (c —>■ 00 ), we have Ep —> 00 and Ep/c^ —> mo so that the 
operator reduces to —moT-i i. This is the non-relativistic operator for the free 
TOA at the origin (c.f. [7J|T2]). Moreover, it can be shown that T is still canon¬ 
ically conjugate to [?7$] = 77$ = Epa^. That is, we still have ['H<s>,T]^ = ih^ 
for an arbitrary $. 
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5 T is a symmetric operator 

For ease of notation, we write T = a^'T ■ If T is symmetric (or Hermitian 
in the sense of m). we should have = ^Td)a|d)f,^ for some 

‘I’a = {4>ai 4’a2)'^ and = {(j)hi (j)b 2 )^■ For simplicity, we let the c^’s vanish at 
infinity so that and are in the Hilbert space Eq. (HU. We first calculate: 

2 2 

pc 

p 

pc^ 

EpC^ 




d 1 

1 d ^ 



^ dp Ep 

= -P^-^Ep = 

Ejdp 

1 

( rriQC^ 

.Ep\ 

1 / pc^rriQC^ 

p dp ' 

K Ep 

+ cE 

~ p[ El 


-TOgc"' + Ep 

p 


2 2 

pc 

p 


Then consider the integral: 

ih 


(j)lT(j)bdp = - 




^ (p4) ) dp 


2„2 


ITIqC 

Er, 


15, d 

-^(Pb + JT 

p op op 



1 

Ep 

1 

p 

(j)b ) dp 




(j)b]]dp 


</>: Ufa dp 


+ 


ih 

T 

ih 

Y 


— oo 
oo 


E„ 


P<^a</’b 

4’bdp + 
{t4>c^ 4’bdp 


Ep^ 

c2 

_5 J_ 

dp Ep 

ih 

Y 


_5 

dp 


P^^ 

_5 

dp 


1_5 

p dp 


(t)bdp 


Er, 


p^c^ p^(? 
P 


£■3 


dp 


dp 


( 28 ) 
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so that we can calculate (noting that Cg = ( 73 ) 




= 


^lasf^bdp 

) 

/ OO 

(t)l2T(j)b2dp 

-00 

(fcj^al) 4>bldp+ [ (T(j)a2) (t)b2dp 

) ' '' J — '' ' 

[[Ual)* (Ua2y 

(crgT^a) ^bdp 
CFs^bdp 


<('61 

<('62 


dp 


= {r<^a 


/ 3 > 

(29) 

where, is just the complex conjugate of the transpose of the matrix A. T is 
then symmetric. This means that its expectation values are real. 


6 The eigenfunctions of T 


Now, we solve for the eigenfunctions of T by considering positive states 
(1 0)^((>+(p) and negative states <b_ = (0 l)^^_(p) so that T<&± = r<i>± 
T is their corresponding eigenvalue. Using Eq (HZl), we get 


$+ = 
where 


T(p± = 


^P'AZ (±<(’±) + ^ (±P<('±) 


2„2 


ITIqC 

2Er, 


2Ep 2 

Er, 


+ 


2c2 


dp 

ih f 1 d 

p dp 


(±<('±) 


dp 


±-<(>± 

P 


±i^^^p'^4'± = +P^4‘± + (2moC^ +/) {^P-^ - 


2 Ejnr) To, q 

±i—^p (()± = 2 p^ 


±i- 


h 

2E„t 


d(l)± 

dp 


2mQC^p 


d(l)± 

dp 


- (p± 


p' 2 I 22 

p + mgC 


d^± 

dp 

d(j)± 

dp 


±i 


PC^T 
Eph 


D± = 


2„4 


2pEl d(l)± 


dp 


ITIqC 


±i- 


+ 




<(>± 


( 30 ) 
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We treat the eigenfunction (j)± as a distribution and solve Eq (15(11) for p > 0 and 
p < 0 (mirroring [19) 1 which gives us 





exp 


h 


&{p\p) 


where, A is some constant and &{p) is the Heaviside step function. We choose 
our eigenfunctions as the sum and difference of these so that our eigenfunctions 
would have definite parity in momentum. That is, 











(31) 


Note that sgn(p) is the sign function. This choice is also made in order to re¬ 
produce the dynamical behaviours of the so-called nodal and non-nodal eigen¬ 
functions of the Confined Time of Arrival operators min] as we will show in 
section ini Explicitly, our eigenfunctions are then 



(32) 

For a given eigenvalue t, 'T then has two degenerate eigenfunctions (p) and 
(p) describing a particle with charge Ae. 


7 Completeness of the 

If the eigenfunctions ^^"^.(p) given by Eq (1551) are complete, then we should 
have the outer product 

/ OO 

^A"r(P)®A"r^(p')t^'r = (To6{p - p') (33) 

A=ini 
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Consider first the integral 


J—cx) y ^p^p' J—oo 

= \A\^^^^^27rhSiEp-Ep,) 

= 27r?i|^P^-r^ [S{p - p') + 6{p + p')] 

tjp \p\c^ 

= — [5{p - p') + 5{p + p')] 

(34) 


and similarly, 

4>^x,lip)4>^\,l* {p')dT = 

} 

= 1^1% 


^ \pp'\c^ 

EpEpi 

' \pp'\c^ 

EpEp! 

i2 Ep [(5(p - p') + S(j) + p')] 


■sgn(p)sgn(p') 


„AiT(Bp-£;p,)/R^^ 


sgn(p)sgn(p')27rM(£;p - Ep>) 


= 27 r;i|A|^—sgn(p)sgn(p')- 


E, 


= 27r/i|A| 


'p 

2b|c E.^ 


p 


Ep \p\c^ 


\p\c^ 

[S{p-p) - S{p + p)] 


= — [S{p-p') - S{p + p')] 

C 


(35) 


where, the second to the last line was obtained since 6 (j) + p') implies that p 
and p' have opposite signs. The left hand side of Eq (l33|) then reads 


E E 


f 0(A) 

l0(-A) 


(0(A) 


/ OO 

-OO 


- /^0(A) 0 \ 

= —\A\^aQS{p-p') 
c 

= cro6{p-p') 




\A\^S{p-p') 


(36) 
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where, 0(A)0(—A) = 0, and we let A = so that the eigenfunctions 





h 

exp ( ) sgn(p) 


( 37 ) 


form a complete set. It can be shown that upon time evolution, the set of 
eigenfunctions Eq dSlD still form a complete set for any later time t. 

Before we proceed, let us consider the non-relativistic limit of Eq (l37l) for 
p > 0 and p < 0. That is, as c —>■ oo, we have Bp ^ mgc^ + and Epjc? ^ mo 
so that 


$ 


(pa) 

A,r 





0(PAP) 


which, apart from a phase factor irrelevant when taking probability densities, 
are just the eigenfunctions of the free non-relativistic TOA operator —moTLi i 

[HiiniiiiiE]. 


8 The Non-Orthogonality of the 

In this section, we wish to compute for the inner product 

/ CX) ^ 

^A"r'(7’)0'3«>y^(p)rfp 

-oo 

= (0(A') J 4'yyip)(t)tl{p)dp 

/ QO 

-OO 

/ oo 

-OO 

(38) 


Note that when A 7 ^ A' (they differ in signs), 0(A')0(A) vanishes so that the 
product (0(A') 0(—A'))(0(A) — 0(—A))^ also vanishes. Also note that when 
n ^ n' the integrand t* ip)4'^xl(p) P integral vanishes. 
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We are then left to compute for the integral in the last equality. 





Anh Ep P 

c 

2TTh 7 q Ep 
1 

It:} 

—— lim 
ZTtR e^0+ 




(r - t') 


Ep dEp 


exp — I e — Ai 


. (r - r') 


Ep j dEp 


I exp moc^^ 


lim 
27Th e^0+ 




1 ^ (T-T')moc2 e + A*—r— 

— hm e^* K - ^^ 

ZTrfi e^0+ 




where exp(—e) was neglected since it is only a multiplicative factor approaching 
unity. From |34] . we recognize a term above as a Dirac Delta. 


4‘^xl*(.p)4>tlip)dp = 


^Xi 


(T-T')r^ 


' ^ A,-(lull 

(5(t -t')-\ -- lim - - - r 


6{t — t') Xie 




2t:(t — t') 


Therefore, Eq. (1381) becomes 

/ OO ^ 

-OO 

— dx^X'Sji^n' X 


5{t — t') Xie 


Xi 


2t:(t — t') 


(39) 


The eigenfunctions are then non-orthogonal. It can also be shown that 

upon time evolution, the set of eigenfunctions remain non-orthogonal at any 
later time t. This suggests that the operator T is non-self-adjoint. T is then 
a maximally symmetric operator. A property which is also true for the non- 
relativistic case — moT-ip [nmniieiiiTiiiiiiii]. 
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9 The eigenfunctions in the ^ — x representation 
and their dynamic behaviour 

Now, we consider the time evolution of Eq (157)) by solving the Klein-Gordon 
equation 

0 t) (40) 

for t) = (0(A) Q{—X))'^(l)^^l{p,t) where is the $ — p representation 

of the Hamiltonian in Eq (fT51) and 0) is given by Eq (1571) . The solution 

of Eq (l40l) is readily seen as 


4’tr (P: t) = </'±(r (p> 0) exp ( ) = ( T 




Er,t 



4>±)r{p^t) = (/'±,r(p>0)exp ^ 


We then write the time-evolved eigenfunctions as 


exp I =Fz 


{t - T)Ep 

h 

{t - T)Ep 


sgn(p) 

(41) 


4>yr(P:i) = (^0(-A) ) 
0(A) 


0(A) ^ H+) 


bk 


0(-A);vi^V^'"'’'"^* 


. {t - T)Ep 




(42) 


We can also rewrite the f)’s by letting q = 


‘fK(9.0 = {if%) - p)=e”(9) 


(43) 
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We then take its Fourier Transform to obtain the eigenfunctions in the position 
representation. That is, 



7S 


Explicitly, we have 


_ / 0(A) \ ^ 

-W-X)) 2^hy2j_^0T7 


0(A) \ rnpc /c 


©(—A )J irh Y 2 
0(A) A moc fc 


</TT< 


: COS 


0(-A); Trh\ 2 Jo ^ 


/TOOC 

V h 

/TOOC 


xq) 


\ 

f moc^ 

xqj 

cos —— 
\ n 




(< - t)\/1 + dq 


— iX 


f 0(A) A WqC /£ 
1,0(-A)j 7r?i V 2 




0 

0(A) A w U(A+) 
e(-X)J nh\l2^-'^’^ 


/moc 
COS I —^xq 
V 


') sin (^ 


moc 


- *A/2^(t^(x,f)) 


(t - t)\/1 + ) dq 


(44) 


^i.T J exp (g, t)c^<? 


\/2TTh J- 

©W ) ^ /I r 

0(-A); 27r/i V 2 y_oo 

e(A) 'i /| r p-A.=#(-)v^d, 


0(—A)y 7r?i 

0(A) A imoc fc 


0(—A )J ttH V 2 
! 0(A) ^ imoc fc 


\Q{—X)J Trh 

0(A) A imo. 
Q(~A)y -kH 


“ y/q 


sin 


moc 




(. 

/TOqC 

V h 


xq cos 


moc 


(t - t)\/1 + g2 dq 


: Sin 


xq^ sin 


moC 




(f-T)\/l + g 2 ) dq 


(45) 


where the /’s are just the integrals and note that they are independent of the 
charge signs A. We can then calculate the charge density in configuration space 
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by using the $ — a: version of Eq. (US 

^ Q(-^)) {\fi?hx,t)\^ + |A/^"^(a;,<)p) 

= (l/}”r + l/2"Ha;,0P) (46) 

which is clearly positive or negative definite. We can then interpret 

^ + \f^:^{x,t)\^) (47) 

as the time evolution of a probability density (in configuration space). Note 
that PT"'\x,t) is the same for A = ±1. We now have the means to see how 
the probability density of a positively or negatively charged particle described 
by the state t) evolves through time. However, the integrals / strictly 

diverge so that in order to calculate Pr'^\x,t), we need to insert a converging 
factor exp (—eg) and let e —>■ 0. In the contour plots, we set h = c = mo = 1 and 
choose e = 0.3 to be sufficiently small. 


j=<+) j=<-) 



Position X Position x 


Figure 1: Time Evolution of the Probability Density with t = 0.5 

As time approaches the eigenvalue r, it can be observed in the pl+1 contour 
plots that the probability density peaks about the origin, while in the 
contour plots, the two probability peaks approach a; = 0, albeit the probabil¬ 
ity always vanishing at exactly a; = 0. Because of these behaviours, we call 
(a;, t) the non-nodal eigenfunctions and (a;, t) the nodal eigenfunctions, 
in accordance to nmi]. The localization observed near t = t and x = 0 can 
be interpreted as the particle being very much likely to be found at the origin at 
the time of the eigenvalue. Thus, if we are to accept that the measurement or 
appearance of a particle at a certain location be characterized by the localiza¬ 
tion of its wavefunction at that location, then the plots suggest that the particle 
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^+) 



Position X Position x 


Figure 2: Time Evolution of the Probability Density with t = 1.0 


being described by unitarily arrives at the position of the origin at 

the time of the eigenvalue much like in the works of Galapon et al m- That 
is, the particle is in a state of definite arrival time. 


10 The TOA probability distribution n 0 ^(r) 

We now calculate the time of arrival probability distribution for a particle de¬ 
scribed by the state = (0(A) That is, we wish to find the 

probability (density) that the state of the particle will ’collapse’ into a state 
of definite arrival time = (0(A) 0(—A))^())^"^. Explicitly in momentum 
representation, 11^^ (r) would be given by 






^\ip)4>^\)r{p)dp 


^\{P) 


-t- 


\/ 47r?i'' 

IIpIc 

1 Ep 


$xip)\^ 

Airh 

\p\Cx^rE 

V Ep 


(48) 


We then essentially need to evaluate the overlap integral between the initial 
state and a TOA eigenfunction Note that from section |51 states with 
X X' have vanishing overlaps. This means that a positively (negatively) 
charged particle state has a vanishing probability to collapse into a negatively 
(positively) charged particle state. 
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As an example, we consider a particle in an initial state given by $ = (1 
where 


Hp) = 




exp 


{p - Pof _ .Xop\ 


which can be shown has expectation values of momentum and position < p >= 
= f-oo ^'•'o' 30 'oP^dp = Po and < a; >= ($ |#)^ = 

Xo and spreads in momentum and position Ap = y/< p'^ > — < p >'^ = and 
Ax = y/< x'^ > — < X >‘^ = , respectively. The probability that the particle 

described by will arrive at the origin during the times r and t + c?t is then 
given by Ii^{T)dT where 


n0(T) = 


+ 


yJmQCy/l^ 

° 1 


{p—p^r I .• ^op 
^2^2 h 

e ^ 0 ^ 


47Th\ E, ® 


Cp— Po) i.-pop 
„2-2 +®“!r" 

g PigC 




(49) 


We set h = c = Too = 1 to illustrate the TO A at the origin probability distri¬ 
bution for a particle with different (average) initial momenta and positions po’s 
and xo’s, respectively. Specifically, we start with a particle moving to the right 
from the left side of the origin. 


Initial Position = -7.0 



Figure 3: TOA Probability Distribution for xo = —7.0 and different po’s 

Consider first the cases where the initial average position of the particle 
Xo is sufficiently far from A^, = One can show that for initial aver¬ 
age momenta po sufficiently far from Ap = the peaks of n 0 (r) corre¬ 

spond to r’s which resemble the free classical relativistic TOA at the origin 
tciass = —xoa/po -b j(poc). That is, the most probable times of arrival at 
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Initial Position = -1.0 


Initial Monnentum = 0.1 
Initial Monnentum = 0.5 
Initial Momentum = 1.0 
Initial Momentum = 1.5 
Initial Momentum = 2.0 
Initial Momentum = 5.0 



Time of Arrival 


Figure 4: TOA Probability Distribution for xp = —1.0 and different pp’s 

the origin r^p’s reduce to the tcZass’s. Note also that for these cases, the r^p’s 
approach (for increasing po) but are always greater than the TOA at the origin 
for a photon tph = —x^jc denoted by the vertical dotted lines in the figures. 
The time tph then serves as a limit as pq becomes very large, as expected for 
a relativistic particle. For being a quantum particle as well, the TOA distribu¬ 
tion also have spreads from their respective maxima which may pass through 
the vertical dotted lines and drop very close to zero as |r| becomes very much 
less than tph- However, this suggests that there is still a non-zero probability 
for the particle to be superluminal. It can also be seen that the peaks become 
higher and sharper as po becomes farther from Ap suggesting that the particle 
becomes more ’classical’. Now consider the cases where po is sufficiently close 
Ap but with xp still far from A^,. One can observe in these cases that the peaks 
have become lower, broader, more flattened, and the Tmp’s become farther from 
tciass suggesting that the quantum behaviour of the particle has become more 
apparent. Specifically, the probabilities become more evenly distributed, slowly 
tapering off to larger r’s, so that the probabilities of the larger times of arrival 
would become just as probable as Tmp for pp ^ 0, still, in a way, supporting the 
classical notion that a slow moving particle would take a longer time to reach its 
destination. In contrast, the probabilities more sharply approach zero for |r|’s 
less than tph, suggesting that superluminality of a particle is still less likely. 
Also, since the momentum probability distribution of the particle can now have 
a significant spread into momenta with directions opposite to po, the state of 
the particle can be interpreted as having a rather significant non-zero probabil¬ 
ity of moving in the opposite direction, reflected in the small probabilities for 
negative r’s. That is, there is a non-zero probability that the particle has al¬ 
ready arrived at the origin. For the case where xp is sufficiently close to A^,, the 
probabilities spread through tph and into the negative r’s even for pp far from 
Ap. This does not necessarily suggest superluminality, however. This can be 
taken into account by the position probability distribution having a significant 
spread at and through the arrival point (the origin), suggesting that the particle 
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has significant probabilities to be initially located to the left, right and even at 
the origin. Also, when po ^ 0 and xg ^ 0, the probabilities become more spread 
out and \Tmp\ become less than tph- The quantum behaviour is most apparent 
because of the significant probabilities of the particle moving to the left or right 
and the particle initially located at the left or right of the origin. 


11 Conclusion 

In this study, we have derived a maximally symmetric, one-particle TOA oper¬ 
ator T canonically conjugate with the Hamiltonian of a free spin-0 relativistic 
charged particle in one spatial dimension. In the non-relativistic limit, we see 
that (for positive states of the KGE) T reduces to —rngT-i i, the TOA operator 
for a non-relativistic particle with mass mg. Being a one-particle operator, it 
does not mix the positive and negative states of the KGE. That is, if we start 
with a particle, we would still end up with the same one. It then makes sense 
to consider systems starting with a relativistic particle in some initial state $ 0 ; 

T$o ) which, as 

/ 3 > 

we have seen, is a real number (as desired for expectation values). We also 
solved for the eigenfunctions of T (which form a complete and non-orhtogonal 
set and reduces to the appropriate non-relativistic limit apart from an irrele¬ 
vant phase factor), and investigated their dynamical behaviour by considering 
the time evolution of their associated probability densities PT"'\x,t) (which is 
non-negative as expected). Note here that PT^\x,t) is the same whether we 
are considering a positively charged or negatively charged particle. We also saw 
from the plots of PT”'\x,t) that the eigenfunctions become more localized at 
the origin at their corresponding eigenvalues. This suggests that a particle in 
a state is in a state of definite arrival time t at the origin as what we 
should expect for the eigenfunctions of a TOA operator. We also calculated 
for the probability distribution n^^(r). That is, we calculated the probability 
(density) that a particle in some initial state would arrive at the origin at time 
T. As an example, we let the initial state of a particle be normally distributed 
about some various average initial positions and momenta, Xg's and po’s, respec¬ 
tively, and plot their corresponding probability densities. As expected, the plots 
revealed relativistic classical behaviour when Xg and pg are sufficiently far from 
their respective spreads, while quantum wave-like behaviour became apparent 
when either Xg or pg are close to their respective spreads. 

These results give us confidence that we are on the correct path in promoting 
time to an observable in relativistic quantum mechanics. This may be in a step 
forward towards having a formalism where physical observables are on equal 
footing, perhaps paving the way in reconciling the inconsistent notions of time 
in quantum mechanics and general relativity. As for some speculative examples, 
upon decomposing or foliating spacetime into spacelike hypersurfaces and letting 
them evolve through a certain time variable which leads to the Wheeler-DeWitt 
equations [35) . perhaps one can define appropriate ’momentum and energy’ 


and asking for, say, its average TOA at the the origin (<i)o 
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operators for the gravitational field and impose that they should be canonically 
conjugate to some ’space and time’ operators whose eigenvalues may be the 
spacetime points. Although interpretations of such operators should be studied 
further. Elsewhere, in Loop Quantum Gravity, the discreteness of space was 
implied through the volume operator |36] . Similarly, perhaps a time operator, 
conjugate to a certain Hamiltonian, can also be constructed which may imply 
the discreteness of time. Or perhaps there is a spacetime volume operator so 
that there may be a smallest chunk of spacetime. Also, in the continuum 1+1 
dimensional model of Causal Dynamical Triangulations, one can calculate a 
Hamiltonian governing the (proper) time evolution of a universe of some initial 
size or length [SHIMIISS]- As the spatial length of the universe is considered as 
an observable in the formalism, perhaps the temporal length of the universe may 
also be considered as an observable which may be represented by an operator 
conjugate to the said Hamiltonian. 

Future studies may include the interacting case. Constructing one-particle 
operators may be difficult, however. This is because separating the equations 
of motion for positive and negative states may not be possible. In the picture 
provided by QFT, the interacting fields suggest that the particles associated with 
these fields collide and interact - creating and annihilating particles and anti¬ 
particles. That is, the non-conservation of particle number may be necessary. 
Thus, restricting our picture to one particle would be more difficult. Along the 
lines of the interacting case, a similar system may be formulated in a curved 
spacetime. The formulation of the equations of motion (the KGE in this case) 
should then be able to accomodate this nontrivial curvature. Another area 
would be to relate possible results here to QFT since they should be able to 
describe the same systems. As an example, both approaches should be able to 
provide a probability amplitude for a particle to start from an initial spacetime 
point and arrive to a final spacetime point. Also, since the derivations in this 
study is formal, further studies on the domain of 'T is needed. We should be 
able to define a dense domain on which “T can act upon for it to be a meaningful 
quantum operator. Lastly, this study focuses on the spin-0 particle. We should 
also be able to construct a one-particle TOA operator for particles with other 
spins (i.e. spin-1/2) canonically conjugate with its Hamiltonian. 


A Some useful relations 

Consider the commutator between the momentum and position operators p and 
(j, respectively, [p, g] = —ih. We calculate the following: 
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[p,9l = [P,?"-'?] 

= [P, q] + [P, q’'~^] q 

= -lhq'^-^+[p,q'^-^]q 
= -ihq’^-^ + [p, % + [p, q^-"^] q)q 

= -mq>^-^+[p,q>^-^] e 


= -ihkq’^-^ 

=p[p,q^] + [p,q^]p 

= -ihkpq^-^ + [p, q^] p 

= -ihk{q’^-^p - ih{k - l)q’'-^) + [p, q'"] p 
= -2ihkq’^-^p - h^k{k - l)q’^-^ 

= p[p,q'^]-lhkq'^-^p 

= P [p, q'"] - ihk{pq'"~^ + iHk - l)f-^) 

= -2ihkpq'^-^ + h^k{k - l)q^-^ 

(50) 


Now, consider the complete and linearly independent set of Bender-Dunne oper¬ 
ators denoted by Tm,n- These are the Weyl ordered quantization of the classical 
monomial p^q^ 


^ OO 

fm,n = -^ X! 

- OO 

= -F 

On / ^ 


n\ 


2" k\{n — k)\ 


n\ 


2" k\{n — ky. 


qkpm^n-k 


^u-k^m^k 


(51) 
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and calculate the following: 


^Tm,n - X! 


'-2 ^n—k 




= -V 

On / ^ 


n! 


{q^p^ - 2ihkf-^p - n^k{k - 


2^ k\{n — k)\ 


— Tm+2,n ^^^^m+l,n —1 


h^n{n — 1) 


2^m,n—2 


T’m.np" - 2„ E ^!(„ _ A:)! 


n\ 


^ oo 

= -y 

on Z—^ 


n! 


2^ k\{n — k)\ 


q'^-kpTq'^f 


q^-kp"^{p^qk + 2ihkpf-k - h^k{k - l)r~ ) 


— Tm+2,n “t” 1 


k?n{n — 1) 


—2 


(52) 


Consider also the linearly independent, complete set of 2x2 matrices, denoted 
by 



Note that cro is just the identity matrix. Using the property of the last three 
matrices UxOy = iaz where, (x,y,z) are just cyclic permutations of (1,2,3), 
and (Jjaj = ctq for j = 0,1, 2,3, we arrive at some relevant commutators and 
anti-commutators of the matrices aj. Namely, 


[(73, CTj] = 2i5j^i(j2 — 2idj^20'i 
[(73 -|- ia'2, cTj] = 2Sjp(cr3 + i<J2) — 2i6j^20'i — 26j^3ai 
{(73 -|- 1(72, CTj} = ‘26jp(a3 -|- 1(72) + 2,i6j^2O'0 + 2(jj,3(7o 

(53) 


Lastly, consider the transform U and its inverse U ^ 
tion. 

jj±i _ ( wqc ^ -H Ep)ao =F ( toqc ^ - 
^JAmoc^Ep 


in momentum representa- 


Ep)ai 


(54) 
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where, Ep = + m^c'^ and calculate 


d ^ pc^ 

dp ^ Ep 


^JJ-I ^ TOoC^(go + CTi) l9 (p-l/2\ (gp - CTi) d 


dp 


\/4moc‘^ dp 


{Ep'^) 


Vdmoc^ dp 




^ moC^{at+ai) f_}_j^-3/2P!_\ " ^i) (}_p-i/2^ 


\/4moc^ 


2^^ E„ 


\/4moc^ 


E„ 


_ I nioc^jai + q-q) _ (cti - cro)Ep ^ f_P^\ 

y ^y4moc^Ep yj4moc^Ep J V ^E^J 

= u-d-fsL) 


CTl 


2i?pV 


CTl 


Ua^U-^ = 


CT3(tooc^ + + 2ia2imlc‘^ - Ep) + (j3(moc^ - Ep) 


4moc^Ep 


2mnC^+p^c^ p^c^ 

= ^3— -tt;-*ct2: 


2moc^Ep 


' 2moc^Ep 


moc" . . P 

= (^3— -h ((73 - 1(72) 




2moEp 


U{a3+ia2)U =Ua3U + 


_i ia2{moc^ + Ep)^ + 2ct3(—+ ia2{moc^ — Ep) 


= Ua^U ^ - (73 


2moEf 


+ 1(72 


4moc^Ep 
2mQC^ + p^c^ 
2moc^Ep 


= (as + 1(72) 


mpc 

E„ 


(55) 


We can then calculate the Klein-Gordon Hamiltonian in the 4>—p representation 


P^ 


= U [ (as + i<y2)-z -1- asmoc U $ 


= U(as + ia 2 )U 


2mo 

-1 


^ + UasU ^TOoc^<i> 


2mo 


, . .moc^ p^ ^ 

= (as + 1(72)^^—::-$ 


Ep 2mo 


moE ■ \ P 

(73 — -h (as - 1(72 ) 


Er, 


2moEp 


moc^4> 


p^c^ + rnlc^ 

= as- - 4> 


Er, 


= asEp^ 


(56) 


Note that the Hamiltonian is even (i.e. it does not mix positive and negative 
states) 
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B Eigenfunctions in the ^ — a; representation 

Using Eq. 63]), we transform it back to the Schrodinger (^k) representation by 

where, 


[7-1 = 


(mpc^ + Ep)aQ + (wqc^ - Ep)ai 
(1 + \/l + q^)(7o + (1 - v^l + q^)ai 

2 yrT 7 


In the — p representation, the eigenfunctions then take the form 


(57) 






\/M 




A + -^x /l + g" 


4\/l + ^2“ Vl - AAi + 

c \/M Az^^^(i-T)\/i+v Al + + A 


nh 4 


c vigl ^-Az:r^(t-T)^i+g2 A + A Al + 


71"fi. 4^/1 + ^2 


1 - Xy^l + q‘ 


sgn(q) 


\/M - 


Az- 


nh 4 




(58) 


To obtain the eigenfunctions in position representation, we take its Fourier 
Transform. Namely, 


1 

4/g(cr,t) = ^J 


moc 


V^TtH J-c 




(59) 


Explicitly, the eigenfunctions in — a: representation are 


25/2^^ 


moc 


, 3/2 


2^/^7vh 


«/ —oo 

I cos(^xg)v^, 


)a/i+ 9^ 'i'^) + A 

(l + g2)-i/2_A 


dq 


(l+g2)-i_ J 


25/27rft y_, 




. fmoc \ ^ 

'Jo ^“ 1 — 


2^/^'Kh 


Xi —ft— 

e ^ 


(*-1-)^!+^ ({^ + q^) 2 + A 
(1 + g^)-i - A 


dq 


(60) 
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